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Abstract 


Neutrinoless double beta decay allows to constrain lepton number violating extensions of the 
standard model. If neutrinos are Majorana particles, the mass mechanism will always contribute 
to the decay rate, however, it is not a priori guaranteed to be the dominant contribution in all 
models. Here, we discuss whether the mass mechanism dominates or not from the theory point 
of view. We classify all possible (scalar-mediated) short-range contributions to the decay rate 
according to the loop level, at which the corresponding models will generate Majorana neutrino 
masses, and discuss the expected relative size of the different contributions to the decay rate in 
each class. Our discussion is general for models based on the SM group but does not cover models 
with an extended gauge. We also work out the phenomenology of one concrete 2-loop model in 
which both, mass mechanism and short-range diagram, might lead to competitive contributions, 
in some detail. 
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I. INTRODUCTION 


Experimental limits on half-lives of neutrinoless double beta decay (0z//9/3) give stringent 
bounds on many Lepton Number Violating (LNV) extensions of the Standard Model (SIVE; 
for a recent review see, for example [l|. Recent experimental results give limits for ™Ge [ 2 I 


and ^^®Xe |3|-l5[ in excess of 10^® ys, which place an upper limit on the effective Majorana 
mass (m^) ^ of the order of roughly (m^) < (0.2 — 0.4) eV, depending on calculations of 
nuclear matrix element [ol-lsl. 

However, from the theoretical point of view it is not a priori clear, whether the mass 
mechanism gives indeed the dominant contribution to the double beta decay rate, and many 



1C 

1 and also in R- 


11, 

12 

and bilinear 


parity violating (^p) supersymmetric theories with both trilinear ^p 1 ^^,, 

^p 0, Q terms. Furthermore, leptoquark models E3 and more recently models with 
colour octet scalars 16| or colour sextet diquarks 17H19j| have been discussed. 


^ (m^) is defined as (mi/) = where the sum runs over all light neutrinos. This is equivalent to 

the (e, e) entry of the Majorana neutrino mass matrix in the basis where the charged lepton mass matrix 
is diagonal. 
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Given that there is such a large list of possible lepton number violating models, is it 
possible to determine which contribution to the decay rate is the dominant one? - 

Perhaps, if double beta decay were to be observed in the next round of experiments and 
either KATRIN 20(] or cosmological data 2ll-l23| also find hints for a neutrino mass scale of 
the order of, say, somewhat larger than (P(O.l) eV, one could claim on the basis of minimality 
that the mass mechanism {my) gives (at least) the most important contribution to the total 
decay rate. However, once upper limits on the total neutrino mass {^my) placed from 
cosmology drop below the level of O{0.1) eV, the question becomes exceedingly difficult to 
answer. 

In that case, from the experimental point of view, there remain only a few possibilities 
to make progress. For exarnple, measurements of the angular correlation between the two 
electrons from Oz//?/? 10|,l2J,l2^ or measuring double beta plus decays 2g offer the possibility 


to identify the Lorentz structure (equivalently the chiralities of the emitted electrons) of the 


LNV processes. However, realistically the SuperNEMO proposal 2^ could only accumulate 


the necessary statistics to identify the Lorentz structure, if the half-live of ®^Se is below 10^® 
ys, while there yet exists no experimental pro pos al with a sufficient sensitivity for 0z//?+/£'C' 
decays to make use of the ideas discussed in j26|. 

From the theoretical point of view, as discussed in |27|, [2^, the amplitude of 0i'(3(3 decay 
can be divided into a long-range and a short-range part. Here, short-range means that all 
particles appearing in the diagrams for double-beta decay are heavier than the nuclear Fermi 
scale, i.e. 0(0.1) GeV. Gurrent limits from Oz//?/? decay then correspond to lower limits on 
the effective scale Alnv of lepton number violation. 


Ar 


LNv= >(l-3)TeV, 


( 1 ) 


where Qes is some mean of the couplings appearing in the diagram and ms and mp are the 
masses of the fields that mediate the Oz//?/? process, see the next section for details. This 
mass scale is testable, at least in principle, at the LHG, and the combination of future LHG 
limits (or a possible discovery, to express it in a more optimistic way) and double beta decay 
data might allow to test man y, b ut maybe not all, of the possible short-range diagrams that 


contribute to the decay rate [2^, [30 


In this paper we take a different approach and study the question, whether the mass 
mechanism is dominant or not, from a purely theoretical point of view. As described above, 
current and next generation Oz//?/? decay experiments will test LNV interactions at the TeV 
scale. Such TeV-scale LNV interactions, on the other hand, appear also in the context of 
radiative neutrino mass models. In other words, a new physics (short-range) contribution 
to 0z^/?/9 decay will always also produce a non-zero neutrino mass. In this paper we discuss 
the relation between possible models for short-range contribution to 0z//3/? decay and the 
neutrino mass-generation mechanism. Our study is based on the complete list of possible 
decompositions of the d = 9 (short-ran ge) double beta decay operator given in |3l|. 
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The general decomposition list given in 3l|, l32|] is equivalent, in principle, to defining all 
models which can give a contribution to double beta decay, and the black box theorem. 
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FIG. 1: Schematic explanation of the black box theorem: The theorem guarantees that, once 
Oi^/3/3 decay has been observed, Majorana neutrino masses will be generated, independent of the 
underlying model, at the latest at the four-loop level. By itself, this theorem does not guarantee 
that the mass mechanism is the dominant contribution to Oi^/3/3 decay. 


see below, guarantees that all these models will produce Majorana neutrino masses at or 
below four-loop order. Our approach therefore basically consists in classifying all the possible 
models contributing to 0z//9/3 decay with respect to the loop level at which they will generate 
Majorana neutrino masses. We can then discuss the expected size of the two contributions 
to 0z//9/3: (1) the d = 9 short-range contribution Od=g and (2) the contribution from the 
neutrino mass mechanism (rUj,) which is radiatively induced in the corresponding models 
and conclude model-by-model, which one of the two is expected to dominate. Given that the 


list of (M| is tree-level complete, our discussion is quite general and covers actually models 


of neutrino mass from tree-level models to 4-loop models. 

Before entering into the details of our work, let ns briefly comment on the well-known 
relation between short-range 0z//9/3 contributions and neutrino masses, i.e., the black box 


theorem 
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The theorem guarantees that, once decay has been observed, neutrinos 


are Majorana particles. However, the black box theorem does not demonstrate that the 
mass mechanism dominates 0u(3l3, since it only guarantees neutrino masses at the level of 
four-loop. Obviously, a four-loop diagram, additionally suppressed by m^m^mg/ALNy, can 
produce only tiny neutrino masses, which are many orders of magnitude below of what is 


required to explain oscillation data 3^. Nevertheless, the black box theorem, together with 
the general decomposition of the d = 9 double beta decay operator published in l3ll . defines 
the basic idea of our current paper. Indeed we hnd that all “models” listed in [Ml produce 
neutrino masses at or below the 4-loop order as demanded by the theorem. 

Babu & Leung 


We comment that our work has some overlap with 36|, 1371 and 38 


have written down all SM invariant AL = 2 operators from dimension-5 {d = 5, the well- 
known Weinberg operator 39|) to d = 11 and showed the relation between the effective 
operators and neutrino masses on the basis of black-box like loop diagrams. The authors 


In 


34l | an extension of the black box mechanism with flavour violation has been constructed. 
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of [3a, l37| discuss then possible ultra-violet completions for several example operators and 
give estimates for the scales Alnv, for which those operators can explain current neutrino 


data. The authors of |38| provide a systematic study of these operators, for one-loop and 


two-loop neutrino mass models, and discuss also which of these could possibly be tested at 
the LHC. However, our discussion differs from these papers in that we are mostly interested 
in double beta decay and its relation to neutrino mass. 

We mention also the work of which pursues the link between the short-range contri¬ 
bution to 0i^/3/3 and neutrino masses, but takes a different approach from ours. The authors 
focus on three types of LNV effective interactions which consist only of leptons and Higgs 
doublets, and list the models in which those LNV interactions simultaneously drive both 
the new physics contribution to 0z//5/9 and neutrino masses. The main difference between 


40| and our work is that it is assumed in 40| that new physics resides in the leptonic sector 


only. 

In our classihcation, we also rediscover several models discussed in the literature previ¬ 
ously, like for example leptoquark models 15||, which can give 1-loop neutrino masses as 
discussed in 41| or 2-loop neutrino masses, as in the model of 19| or the one in We 
do not cover, however, the possible contributions from light sterile neutrinos. There exists a 
vast amount of papers on this subject in the literature already 43l-l55| and we have nothing 
new to add on this particular subject. 

Since neutrino mass models must not only produce the correct absolute values of neutrino 
mass, but also reproduce the observed flavour structure of the neutrino mass matrix, one 
also has to pay attention to constraints from flavour physics observables. In the authors 
applied the hypothesis of “minimal flavour violation” (MFV) to effective operators that 
contribute to 0z//3/3 and found that the MFV assumption constrains the effective operators 
to be smaller than the detectable level. In this work, we do not adopt any such theoretical 
assumption on the flavour structures of the parameters in the models. Instead we simply 
consider bounds on lepton flavour violating observables as constraints on model parameters. 
We believe this to be the correct approach since any of the “exotic” contributions to 
decay requires the introduction of new scalars, not present in the SM, with their own Yukawa 
interactions with SM fermions. Thus the whole concept of MFV is not very well funded in 
any of the models of interest for decay. 

A few disclaimers might also be in order here. Our analysis concentrates on the true 
d = 9 operator, i.e. it covers only the short range part of the 0z//9/5 amplitude. Our results 
thus do not cover, for example, the long-range diagrams of R-parity violating SUSY 58| 
or leptoquark models 1^. Also, we limit ourselves to scalar exchange, thus models with a 

are not considered. 
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coupling between new scalars and the SM gauge bosons, such as 
Also, this restrictions implies that we do not cover models with an extended gauge group 
either, especially we do not discuss models with left-right symmetry. And, hnally, the list 
of decompositions in 3]| is complete only at tree-level. Thus, we do not consider cases in 
which the neutrino mass is generated at some higher loop level, while the amplitude 


appears at one-loop order, as for example in the recent papers |M 
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The rest of this paper is organised as follows. At the beginning of section |TT1 


as a 


31 


and discuss some generalities useful 


preparation, we will summarise the main results of 
for the latter parts of the paper. We will then discuss the classification of the different 
possible models and estimate in each case the relative size of the contribution from the 
mass mechanism and the short range part of the amplitude. In section IIIIl we will discuss 
one concrete two-loop model of neutrino mass in more detail. Section IIVI summarises and 
discusses our main hndings. Tables with lists of the different models, classified as described 
in section m are deferred to the appendix. 


II. SETUP AND CLASSIFICATION 

In this section, we classify neutrino mass models based on the decomposition of the d = 9 
0i^/3/3 effective operators, according to the number of loops in the resulting neutrino mass 
diagrams. In each class, we will compare the size of the two contributions to the 0z//9/5: (i) 
d = 9 operator itself, and (ii) the mass mechanism (m^) induced by the d = 9 operator. 
This classihcation therefore allows to identify those models, for which one can expect non¬ 
standard contributions (beyond the ordinary mass mechanism) to be important for On/3/3 
decay. Note, that (m^), is equivalent to the e-e entry in the neutrino mass matrix, (Mj,)ee, 
in the basis where the charged lepton mass matrix is diagonal. Thus, in this section we 
concentrate on the comparison of short-range contributions to the size of this entry in {M^). 
Of course, a complete fit to all neutrino data will need to take into account also all other 
entries in {My). A specihc example, how this can be done and the additional constraints 
from both, oscillation data and lepton flavour violating decays, is discussed in section 3. In 
the discussion in this section, we always keep generation indices in the unknown couplings 
of the different models. Other indices could be constrained combining double beta decay 
with, for example, oscillation data. A specific example for this is worked out in section 3. 


Generalities 


The short-range double beta decay operator 


Od =9 oc uu dd ee 


( 2 ) 


can be generated at tree-level via only two topologies shown in Fig. [2j The bosons (depicted 
with dashed lines in the diagrams) in these topologies could be either scalars or vectors, but 
we will consider only scalar exchange here. Assigning the outer fermions with either the left 
(“L”) or the right (“i?”) chirality in all possible permutations allows to derive the complete 
list of “decompositions” (or proto-models) that can contribute to the Oz//3/3 decay amplitude 


at tree level 
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The fermion propagator in topology-I contains two terms, 

i) +m^ 

9 2 ’ 

- rn^ 


(3) 
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FIG. 2: The two possible tree-level topologies contributing to the Oz^/3/3 decay rate. The outer 
lines represent the six SM fermions, while for the virtual particles appearing in the inner lines we 
consider only scalar-fermion-scalar (left, topology-I) and scalar-scalar-scalar (right, topology-II). 


but in the short-range part of the amplitude the hrst term is suppressed relative to the 
second by a factor of /m^ ~ pp/rn^, where Pf is the typical Fermi momentum in the 
nucleus and the mass is suppossed to be larger than (9(100) GeV. Considering then only 
decompositions which pick the mass term from the propagator results in a total of 135 possi¬ 
ble decompositions for Topology-I (T-I in the following), while there are 27 decompositions 
in Topology-II (T-II), if we limit ourselves to scalar exchange [3l|, l32|]. For tables showing 
the different decompositions see the appendix. 

Babu & Leung jSGj have listed AL = 2 operators from d = 5 to d = 11. Among the d = 9 
operators in their list, the following hve are relevant for double beta decay: 


Oil — LLQdpQdji, O12 — LLurQurQ, Ou — LLupQQdp, (4) 

Oig = LepQdRURdR, O20 = LcRURQuRdn. 


Here, we have suppressed the indices of generation, SU{3)c, SU{2 )r, and Lorentz spinor, 
which are contracted appropriately. In the list of decompositions shown in 3l| , there appears 
also one operator not given in 3q, which is^ 


O- — Cr Cr URdRURdR. 


(5) 


As the black box theorem demonstrates, one can obtain the Weinberg operator by con¬ 
necting the quark legs in these effective operators with the SM Yukawa interactions. From 
the effective operator point of view it seems that neutrino masses are generated at 2-loop 
level from On, O 12 and Om The operators Oig and 02g need an additional SM Yukawa 
interaction with a charged lepton to generate a neutrino mass term, thus they end up with 
3-loop diagrams. The operator leads to neutrino masses only at the 4-loop level, which 
is equivalent to the original black box diagram. However, the classihcation of the neutrino 


^ As discussed below, this operator induces neutrino mass only at the four-loop level, 
reason it was neglected in |36l| . 


Probably for that 


7 














(Si,2,1/2) (81,2,1/2) 

L/ei) L(eL) ' ' 


I I 



FIG. 3: To the left: Diagram for Oi^/3/3 decay via charged scalar exchange for Babu-Leung operator 
Oil (BL^^ll). To the right: Tree-level neutrino mass generated via seesaw type-I, using the same 
vertices as in the diagram on the left. Here and in all Feynman diagrams below, arrows on fermion 
lines indicate the flow of particle number, not the chirality of the fermion. 


mass models with respect to the number of loops, which we discussed in the introduction, is 
modihed from this naive expectation, once the decomposition of the operators are specihed. 
In fact, as shown below, many decompositions of the operators in Eq. (j4]) contain auto¬ 
matically the particle content (and interactions) such that neutrino masses are generated at 
lower order, i.e., both tree-level and 1-loop neutrino mass models are found. And, surpris¬ 
ingly, also the opposite case exists: If we restrict ourselves to decomposing the operators 
only with scalar and fermion mediators, none of the decompositions of the operator 
generates a genuine neutrino mass diagram at the 2-loop level. We will come back to this 
important point later in more detail. Genuineness is one of the key concepts in our classih- 
cation method. The term, genuine n-loop neutrino mass model, is dehned as the model in 
which the neutrino mass term is generated at the n-loop level and for which, simultaneously. 


diagrams with loop level lower than n are guaranteed to be absent, see 63|] for more details. 


B. Tree-level neutrino mass models 

Let us start with a rather trivial example, illustrated in Fig. [3l in which the d = 9 
contribution to On/3/3 is related to the Majorana neutrino mass at the tree level. Here, as 
everywhere else in this paper, subscripts on helds denote their transformation properties (or 
charge in case of f/(l)y) under the SM group, SU{3)c x SU{2)i x f/(l)y. A new scalar held 
is denoted by the symbol S, and a fermion held, which is understood as either a vector-like 
fermion or Majorana fermion, is denoted as '0- Thus, "01, 1,0 has the same quantum numbers 
as a right-handed neutrino, while 81 ^ 2 , 1/2 is equivalent to a (copy of) the SM Higgs doublet. 

The Lagrangian producing the left diagram of Fig. |3] necessarily contains the following 
terms: 


^ {y^SQd)ij Qi ■ 5 * 1 ^ 2 , T Li ■ ^p(0’l,l,o)^01,l,O T H.C.. ( 6 ) 

Here, the singlet fermion held 01,1,0 is allowed to have a Majorana mass m^. The dot (•) 
denotes the anti-symmetric tensor (ir^) for SU{2)l. This Lagrangian generates an ehective 









operator for a short-range contribution to 0z//9/3, 


C.d=9 — 


{{ySQd)ll{ySy^)ef x -y r c 




(Q^r) • LL'^ • (gci,^) + H.C., 


(7) 


which corresponds to Of^ in the notation of |28l| . (Here we use the notation for the hve 
relevant short-range operators, dehned in |28| , to distinguish them from the lepton number 
violating operators Oj, j = 11,1^,^19,20 and ‘^ 8121/2 nms^ of S'!, 2 , 1 / 2 - 


Following the method adopted in 


Q 
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and using the experimental bound 


j.0j^/3/3(136Xe) > 1.6.10^5 [ys] 

one hnds the bound on the coefficient of the d = 9 operator as 

1 ((hgQd)ll(F5t/V>)e)^ ^ Gp ^{RR}R ^ Gp ^ 

8 ^5121/2’^^ 2 mp ^ 2 mp 

which can be interpreted as a constraint on the Yukawa coupling: 


< 1.5 • 10-6 


1.0 \ / ”^5i.2,1/2 V / \ 

V100 [GeV]J Vl00[GeV]J- 


( 8 ) 

(9) 


( 10 ) 


If the scalar mediator S'!. 2 , 1/2 acquires a vacuum expectation value (vev), the right diagram 
in Fig. 131 which has the same topology as the type-I seesaw mechanism, will contribute to 
neutrino mass(es): 

rriu) = -- 

my 



The experimental bound on the effective neutrino mass (m^) < 0.3 eV, which is found from 
Eq. dH]) under the assumption of the mass mechanism being dominant, gives 


(Xsiyip'^'’ y 1-0 


10 


-12 


A // jrnj) 
100 [GeV]y V(F)y V0.3[eV] 


( 12 ) 


where we have used vsu — 174 GeV. If 5'i^2,i/2 is identihed as the SM Higgs doublet FT, as 
in the ordinary type-I seesaw, the constraint on shown in Eq. flT^ is obviously much 
stronger than Eq. ffTOj) .^ In other words: The mass mechanism dominates the contribution 
to Offd/d decay by far, if we consider SM Higgs exchange. 

However, -S' 1 ^ 2 , 1/2 is not necessarily the SM Higgs, it could be an additional new state, 
such as appear, for example, in multi-Higgs doublet models. In this case, neutrino masses 
would still be generated through the type-I seesaw mechanism, with the vev (S) of the scalar 


^ If S'!, 2 , 1/2 is the SM Higgs doublet, the coupling YsQd is identified as the down quark Yukawa coupling, 
Yd ~ 3 X 10“^, and the constraint shown in Eq. (ITUl) actually becomes less stringent than even the ordinary 
perturbativity bound. 
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FIG. 4: Three different 1-loop diagrams for neutrino mass 78|, which can appear in one of the 
Oi^/3/3 decay decompositions. For a discussion see text. 


•S'!, 2 , 1/2 independent of the SM vev wsm, such as occurs, for example, in the neutrinophilic 
neutrino mass model of 64-67 . For this case one hnds that if the relaxed constraint 


(5) < 0.14 [GeV] 


{YsQdhiV flOO [GeV] 


1.0 


ms 


(13) 


holds, Eq. (ITOll becomes more stringent than Eq. (IT^ . i.e., the short-range diagram will be 
the dominant contribution to the Oi//3/3 decay amplitude in this case. 

If 5 'i, 2 ,i /2 has exactly zero vev, in the literature often called the “inert doublet”, we can 
no longer directly relate the the relative size of the d = 9 operator with the mass mechanism. 
The only conclusion one can derive in this particular case is the trivial constraint that the 
standard model Higgs coupling with L and '0i,i,o niust obey Eq. (HI. 

As is well-known, there are only three types of tree-level mass generation mechanisms 
(seesaw mechanisms) called type-I 68-71], type-II 0 , and type-III [^. These are 

mediated by the singlet Majorana fermion 'ipixo (type-I), the triplet scalar Si ,34 (type-II), 
and the triplet Majorana fermion " 01 ^ 3^0 (type-HI). From the complete list of decompositions 


given m [M|], one can hnd that T-I-l-i, 2-i-b, 2-ii-b, 2-iii-a, 4-i, and 5-i contain the relevant 
fermion mediators, and T-I-l-ii-a, 1-ii-b, 3-ii, 3-iii, T-II-1, and T-II-3 do contain the scalar 
mediator. From the discussion above, we can conclude that for all of these the short-range 
contribution will be much less important than the neutrino mass mechanism, unless the vev 
of the new scalars 81 ^ 2 , 1/2 or S '!,34 are heavily suppressed compared to usm- 


C. 1-loop models 


As shown in 78| (see also 79|), there exist a total of four genuine 1 -loop diagrams, which 
can contribute to the d = 5 Weinberg operator at the renormalizable level.® Three diagrams, 
shown in Fig. 01 can be related to Oujdld decay decompositions.® Here, we have added a ‘V” 


® There are also three more non-genuine diagrams, discussed in 781, which can be understood as one-loop 
generated vertices for one of the three tree-level seesaws. 

® The remaining diagram Tj^-l-i can be understood as opening-up of the quartic scalar vertex in Ti/-3, by 
inserting an additional scalar. 
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to the naming conventions of 78(], in order not to confuse the 1-loop neutrino mass diagrams 
with the double beta decay topologies. 

Let us discuss the relation between the neutrino mass diagram Tz/-l-ii, which is shown 
as the left-most diagram in Fig. 01 and the decomposition of the d = 9 diagram 

(uLdR){'^dii)(uL'^), which is classihed with the ID number T-II-2 On in 


31 


as an 


example of this class of models. The decomposition leads to a Lagrangian, which contains 
the following terms: 


£ = (VL,shU ■ (S3,2,l/6)^dn,j + (VLQshQj ' ^^3,1,-1/3 (14) 

+ 0^Qds)ijQiSl,2,l/2dRJ -|- /^555('S'3,i_i/3)1S'3^2,1/6'S'J 2,1/2 + 

+ ^^g(S'3,2,l/6)^<S'3,2,l/6 + ^ (S'3,i,_i/3)lS'3,i,_i/3. ^ ^^^(S'l,2,l/2)^<S'l,2,l/2• 

The lepton number violation can then be assigned to be due to the presence of the coupling 
fisss- With this Lagrangian, Eq. (1T4)) . the effective d = 9 Lagrangian that contributes to 
Qufdfd process as short-range effects is given as: 


Cd =9 = (LdR) ■ {QdR){Q ■ L") + H.c. 


9 2 2 

m% mt mo 

*^ 3 , 2 , 1/6 *^ 3 , 1 ,- 1/3 *^ 1 , 2 , 1/2 


D 


(XLds)el(XQds)ll(XLQs)le/^SSS 1 f^nSR\ 

){RR}R 


rriq ran ran 

* 3 , 2 , 1/6 * 3 , 1 ,- 1/3 * 1 , 2 , 1/2 


16 


(15) 


The experimental bound, see Eq. ([8]), on Oa/d/d decay can then be interpreted again as an 
upper limit on the new leptonic Yukawa interactions as: 


(XLds)el(XLQs)le < 0.15 


2 2 2 
ran ran ran 

* 3 , 2 , 1/6 * 3 , 1 ,- 1/3 * 1 , 2 , 1/2 


1.0 [TeV® 


T.0[TeV] 


1.0 


(4Q(is)ll 


(16) 


With the interactions shown in Eq. (1T4)) . neutrinos acquire Majorana masses through the 
diagram Tz/-l-ii. This neutrino mass generation mechanism through the leptoquark-Higgs 


coupling was hrst proposed in |80| and discussed in detail in j^, [81 


Assuming the coupling fisss is smaller than the average of the leptoquark quark masses 
^LQ = { 1 RS 3 21/6 + "^*31 - 1 / 3 )/ 2 ) we can roughly estimate the neutrino mass as 


[ra 


v)a0 — 


'l^dkl^SSS 


IGvr^ 


ra 


LQ 


(5'i/2) 0 ^LQs)dk{ylds)ka + {/d -H- a) 


(17) 


where rad^. is the mass of the down-type quark (of generation k), which enters this estimation 
as the vertex at the left-upper corner of Ta-l-ii of Fig. 01 Applying the bound on the effective 
Majorana mass (ra^) < 0.3 eV (which is obtained from Eq. (IH]) with the assumption of the 
mass mechanism dominance) to Eq. (ITTj) . we have 


(XLds)ek(XLQs)ek < 0.034 


/4.0[GeV]\ 

V rad, ) 


( rai^d y/1.0[MeV]\ 
y.OlTeVly* V k^sss ) 


/174[GeV]\ / (m,) \ 

V (A 1 / 2 ) Ao. 3 [eV]y 

(18) 
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Assuming that the flavour structure of the new Yukawa interactions in Eq. ffl^ is not 
strongly hierarchical, one concludes that third generation quarks give the largest contribution 
to the neutrino mass. If the vev («S'i/ 2 ) is as large as the SM Higgs vev, this constraint is more 
than six (three) orders of magnitude more stringent than Eq. flTbl) for fc = 3 (/c = 1). Note 
again that if the vev (*S'i/ 2 ) vanishes then, as in the tree-level case, the d = 9 contribution 
to Oi^/3/3 and the mass mechanism in this class of models are independent of each other. 

If the leptoquark mass is set to 0(1) TeV and Yukawa couplings are taken to be 0(0.1), 
the trilinear coupling fisss must be 0(100) keV to reproduce 0(0.05) eV of neutrino masses, 
which is the minimum value necessary to reproduce data on atmospheric neutrino oscilla¬ 
tions. Such a small value of the coupling ^sss can (obviously) be probed only in a Offd/d 
process dominated by the mass mechanism, since there is no other LNV process, for which 
experiments have even remotely comparable sensitivity. 

To finish this discussion of the one-loop neutrino mass case, recall the possible two-loop 
contributions to neutrino masses in this model. As shown in the next subsection (and the 
list in the appendix), the model described with Eq. (lT4ll generates neutrino masses at the 
two-loop level, even if the value of (5'i/2) is identical with zero. However, the one-loop 
contribution discussed in this subsection can easily (even with a value of (*S'i/ 2 ) smaller than 
MeV) dominate over the two-loop diagrams, as can be seen from Eqs. flTT)) - (fT8|) . Therefore, 
we classify this type of the models separately from the genuine two-loop models, which will 
be defined and explained in detail in the next subsection. 

Very similar arguments can be applied to the other two topologies, Tz/-l-iii and Tz/-3, 
shown in Fig. IH Such one-loop neutrino mass models appear in quite a large number of 
decompositions. We give the complete list of this class of models in Table UTI in the appendix, 
together with the additional interaction that is required (and is allowed by the SM gauge 
symmetries) to generate the corresponding one-loop diagram. 


D. 2-loop models 


The effective operators On, On, and Ou contain two lepton doublets, and thus naively 
one expects them to generate neutrino masses at 2-loop level,^ when their quark legs are 
connected with the SM Yukawa interactions. However, this naive picture has to be modified, 
once the possible decompositions of the effective operators are taken into account. 

In this subsection, we will first demonstrate why some of the decompositions of the 
operators On, Ou, and Ou do not generate neutrino masses genuinely at the 2-loop level. 
Here, we use the terminology genuine 2-loop diagrams, following 63[, to imply that the 


^ When the two lepton doublets are anti-symmetric in SU(2)l, an additional loop is necessary to obtain a 
neutrino mass term, i.e., the resulting neutrino mass diagram contains three loops, which correspond to 


Oiia, Oi 2 b and Oua in the list of (37|. In Ovflfd decay only the operators Onb, Oua and Oub, which 
contain symmetric pieces in SU(2)l, can appear. These are called On, Ou and Ou for brevity here and 


3ii. 
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FIG. 5: Decompositions of the Babu-Leung operator 7^14, O 14 oc L L Qdpiu^Q, with two scalar 
{S and S') and one fermion ('tp) mediators. For a neutrino mass diagram the quark legs ur and 
Q must be connected via the SM Yukawa interaction The resulting loops (shown on the right) 
are infinite one-loop corrections to the scalar mass terms, that therefore must be contained in the 
tree-level Lagrangians of the corresponding models. 

corresponding model (or decomposition in our case) generates neutrino masses at 2-loop 
order and that no lower order diagram exists. 

Let us first discuss the decomposition of Ou. This operator contains: 

Lai Lq, Ur , Qa^ Qa'> d,R . (19) 

Here we explicitly wrote the 2-component spinor indices for the fermions: dotted for a right- 
handed held (complex conjugate of left-handed held), undotted for a left-handed held. As 
everywhere else in this paper, we restrict the discussion to fermions and scalars as mediators. 
For Oi 4 this implies that the spinor indices on Q and mr must be contracted for the ehective 
operator being a Lorentz scalar. There are only two choices to assign these two quarks to 
the outer legs of a d = 9 tree diagram: (i) They form a Yukawa interaction with a scalar 
mediator, i.e., (urQ ■ S). This is shown as the upper diagram of Fig. [5l Or: (ii) Each of 
these quarks forms a Yukawa interaction with a fermion mediator -0 and one of the scalar 
mediators, i.e., (llR'ijjS){'ipQS'). This is shown as the lower diagram of Fig. [5l When the 
loops are closed, as shown on the right of Fig. [5l via the SM Yukawa interaction yu{Q-H*UR), 
the resulting quark loop is divergent. In other words, these loops are inhnite corrections 
for (i) a mass term mixing the scalar mediator and the SM Higgs doublet or (ii) 

a corresponding term mixing two scalar mediators m‘gg,SS'. Therefore, the original tree- 
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FIG. 6: Examples of non-genuine 2-loop neutrino mass diagrams, based on the decompositions of 
T-I-l-i Oi 2 (left) and T-I-2-ii-a On (right). For a discussion see text. 


level Lagrangian generating these operators mnst contain these scalar mass terms as connter 
terms for the inhnities, and the qnark loop appearing in nentrino mass diagrams mnst be 
snbstitnted with those scalar mass terms. 

This simple argnment actnally holds for any decomposition of Ou- In short, one cannot 
constrnct a genuine (irredncible) 2-loop diagram for nentrino masses from the decomposition 
of Oin if fhe mediators are restricted to scalars and fermions. For this reason, in the 
appendix the decompositions of On are instead either listed nnder “tree-level” or “l-loop”, 
depending on the additional interaction necessary. 

The next qnestion we mnst address is: Can all the decompositions of On(b) and Oi 2 (a) 
genninely generate nentrino masses at 2-loop level? The answer is no, and the argnment for 
those cases is very similar to the one presented above for Ou. Two concrete examples are 
shown in Fig. |6l Consider hrst the diagram on the left, based on decomposition T-I-l-i On- 
The inner loop in these classes of nentrino mass diagrams actnally corresponds to an inhnite 
(1-loop) correction to the scalar qnartic interaction Xss'hhSS'HH. In the diagram on the 
right, the “inner” loop involving S generates a Ynkawa interaction yq-tpH'^Q ■ H. Again, this 
correction is inhnite, thns reqniring a tree-level connter term which mnst be contained in the 
original Lagrangian. Given these additional (bnt reqnired) interactions, the models contain 
nentrino mass diagrams at 1-loop order. Note that, the left diagram in Fig. [6] corresponds to 
Diagram of Fig. 14 in 38|] (with appropriate Hi ggs insertions) and also Diagram (c) of 

Qnite a nnmber of 
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The right diagram is Diagram (B) in 38| and (d) in 
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Fig. 5 in 

possible decompositions of On and O 12 follow this pattern and are thns listed in the tables 
in the appendix as 1-loop models, together with the additional-bnt-necessary interactions. 

After hltering ont all decompositions that resnlt in non-gennine 2-loop nentrino mass 
diagrams, we have fonnd that for all remaining decompositions there are only three types of 
gennine 2-loop diagrams, all of them based on On- These are shown in Fig. [3 The naming 
scheme in this hgnre follows 63|. In the appendix, we present the complete list of the 


gennine 2-loop nentrino mass models and specify the class of nentrino mass diagrams, into 
which each model falls. In Table IIVI two of the decompositions based on the BL operator 
OiQ are also listed. These appear in the table for two-loop models dne to the fact that the 
intermediate fermion is of Majorana type, i.e., for these decompositions, the “asymmetric” 
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FIG. 7; Genuine 2-loop neutrino mass diagrams based on decomposition of the Babu-Leung op¬ 
erator On- From left to right the diagrams are identified as CLBZ-1, PTBM-1 and PTBM-4, 
following the classification of [^ . The left diagram (GLBZ-1) was discussed in the context of a 


neutrino mass model in [1^, and corresponds to the diagram in Fig. 10 in 3^ and (e) in 82l ]. 


The diagram (PTBM-1) in the middle was discussed in 8^, and corresponds to Diagram (D2) in 


14 in 3^ and (f) in Fig. 5 in 8^. The diagram on the right corresponds to Diagram (G) in 
38 1 and (f) in Fig. 5 [^. A model based on this diagram will be discussed in Sec. IIIII 


operator Oig oc Le^ is always accompanied by the “symmetric” operator On oc LL, and the 
associated operator generates nentrino masses at the two-loop level. The catalogue of the 


effective operators appearing with their associated operators is given in the tables of M 


Among these three genuine diagrams, neutrino mass models based on the CLBZ-1 and 
the PTBM-1 diagrams have already been studied in the context of the decomposition of 


the d = 9 operators 
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Therefore, we will discuss a 2-loop neutrino mass model that 


is associated with the remaining possibility, i.e., the PTBM-4 diagram. Here, as in the 
previous subsections, we compare the d = 9 contribution to with the mass mecha¬ 

nism contribution and postpone the detailed discussion on phenomenology of this model till 

Sec. [ml 

The example we choose is based on T-I-4-ii-b, On- The Lagrangian for this decomposition 
contains the terms 


+ (Xipds)i{'>p6,2,l/edfi^i) 83^2,1/6 + {yLds)ai{Ladji^i) ■ 5'3 2 ^/g. 


( 20 ) 


We use the notation S' 6 , 3 , 1/3 = (*S' 6 , 3 ,i/:^(Tg)fj and V'egp/e = {'^ 6 , 2 ,i/ 6 )x{Te)fj- The tensors 

Here, r is the Pauli matrix vector for a triplet of 


31 


Tg and Tg in the SU{3)c are given in 
SU (2) L. The effective d = 9 operator resulting from this Lagrangian can be written with the 
following linear combination of the basis operators of the short-range contributions 

to 012 (3(3 decay as: 


= - ^^QQs)^^^'^Lr^s)eiY^ds)l{YLds)el . gc) . (2,/rf^)] + R.C. 

_ 


m 


6.3,1/3 


’ 3 . 2 , 1/6 


D 


(hQQs)ll(hLi/.s)e(hpds)l(LLd5)el 


ml 


^ 6 , 3 , 1/3 


mt 


53,2,1/6 




— {0^^){RR}R - —{02^){RR}R 


( 21 ) 
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and the experimental bonnd Eq. ([8]) constrains a combination of the coefficients 


(XQQs)nO^Ltl}S)e{yi}ds)l(XLds)el < 6.3 • 10 ^ 


rric 


’ 3 , 2 , 1/6 ‘’ 6 , 3 , 1/3 

1.0 [TeV® 




(22) 


On the other hand, the nentrino mass generated from the 2-loop diagram based on the 
effective operator Eq. fl2T]) also contributes to 01 /13(3 through the mass mechanism. The size 


of 2-loop neutrino mass diagram can be roughly estimated as 36l-l38l. l63 




Nr. 


mt 


(167r2)2 Alnv 


[{YQQs)33{YL^s)a{yipds)30^Lds)/33 + {(3 ^ a)] 


(23) 


Applying the experimental bound from 0i/(3(3 to Eq. (]23|) and substituting Nc = 6 because 
of the colour sextet combination in the loop, we can place the bound on the couplings of the 
third generation quarks: 


{YQQs)33{YLips)e{Yi)ds)3{YLds)e3 < 7.2 ■ 10 ® 


4.18[GeV] 


rUb 


2 



As this rough estimation shows, the mass mechanism and the short-range part of the am¬ 
plitude give similar contributions to 0i/(3(3 decay. Note that assuming flavour democratic 
Yukawa couplings, the constraints Eq. fl2^ and Eq. fl2T|) become equally strong if the mass 
scale of the new particles is taken to be roughly ~ 300 GeV. For larger mass values, the short 
range contribution can dominate only if Yukawas with index “3” are smaller than those with 
index “1”, otherwise the mass mechanism dominates. A more detailed discussion using the 
full expression for the two-loop neutrino mass integral will be presented in Sec. IIIIl 


E. 3-loop models 

From the point-of-view of effective operators, the Babu-Leung operators Oig and O 20 
require three SM Yukawa interactions to generate neutrino masses: two quark Yukawa inter¬ 
actions and one charged-lepton Yukawa interaction, to convert cr in the effective operators 
to L for a neutrino mass. This fact leads us to three-loop neutrino mass models. However, 
some of the possible decompositions of Oig and O 20 contain the ingredients to generate neu¬ 
trino masses at a level lower than three-loop. In such a case, the lower loop contributions 
can easily dominate neutrino masses and make the contribution from a three-loop diagram 
sub-dominant. This can happen for two reasons. First, there are decompositions based on 
Oig or O 20 , in which On necessarily also appears. We call this “associated” operators and 
classify those decompositions in the class corresponding to those lower loop levels, which 
usually will dominate over the 3-loop contribution. Decompositions T-I-2-iii-a and T-I-5-i 
of Oig with the Majorana fermion '08,1,0; which are listed in the table of two-loop models 
(Tab. HV)) . are categorised in this class. And, second, there are decompositions for Oig or 
(920; in which the 3-loop diagrams are not genuine in the sense that one of the sub-diagrams 
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Q(dL) e_R 



FIG. 8: The Feynman diagrams for the decompositions T-I-4-ii-a (left) and T-I-4-ii-b (right) for 
BL operator O 2 o- The former leads to 2-loop d = 7 neutrino masses, while the latter is an example 
of a 3-loop neutrino mass model, see text. 

I'L ’'L 




FIG. 9: Examples of 3-loop diagrams for the decompositions T-I-4-ii-a (left) and T-I-4-ii-b (right) 
for BL operator O 20 , see Fig. [HI The diagram to the right corresponds to a genuine 3-loop model, 
while the one to the left is not. See text. 

corresponds to the 1-loop generation of a certain vertex. We will discuss this case in a bit 
more detail. 

Take the examples of the decompositions T-I-4-ii-a and T-I-4-ii-b, both O 2 o- The Feyn¬ 
man diagrams are given in Fig. [HI and show T-I-4-ii-a —?■ (mr MR)((iL)(eR)(erd/j) and T-I-4-ii- 
b —t (mr MR)(eR)(dL)(eTdi?) graphically. As we will see, despite the similarity between these 
two cases, T-I-4-ii-b will lead to a genuine 3-loop model, while T-I-4-ii-a will not. Consider 
the examples of 3-loop diagrams for these two decompositions shown in Fig. [9l First of 
all, note that the loop diagrams shown in Fig. [9] should be understood as examples only, 
because there might be more than one diagram contributing to the full neutrino mass matrix 
for each decomposition. In the diagram for T-I-4-ii-a (left), one sees that the innermost loop 
effectively generates the vertex MR'03,2,7/6-f^^ at 1-loop order. This one-loop sub-diagram 
is inhnite and, therefore, a tree-level counter term is necessarily to be included in the La- 
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FIG. 10: Example of a 2-loop d = 7 diagram for the decomposition T-I-4-ii-a. This diagram 
corresponds to Diagram (a) in Fig. fT51 

grangian to absorb the infinity. In fact, the qnantum numbers of the particles involved in 
the loop are such that the necessary vertex actually cannot be forbidden at tree-level by 
the SM gauge symmetry. This tree-level coupling has a value that is not fixed by the Qufdfd 
decay amplitude, but the 2-loop (d = 7) diagram that results from this coupling, see Fig. [TOl 
can easily dominate over the 3-loop diagram.® We have classihed therefore all these cases 
of Oig and O 20 , where such a tree-level vertex is allowed, as “2-loop d = 7” (Tab. |V]) in 
Appendix. The loop diagram based on the decomposition T-I-4-ii-b, which is shown as the 
right diagram of Figure [9l does not contain such an inner loop and, thus, such a construction 
is not possible for this decomposition. Decompositions of this type are therefore classihed 
as genuine 3-loop models in the appendix. 

We give here a rough estimate of the size of the neutrino mass generated by the 3-loop 
diagram based on T-I-4-ii-b in order to present some general arguments on the relative size 
of the d = 9 contributions and the mass mechanism contributions to Oufdfd in this class of 
models. This example, in which the d = 9 is mediated by a di-quark, S' 6 q, 4 / 3 , a leptoquark, 
*S'3,2,l/6, and an exotic colour-sextet fermion, 'ipe, 1 , 1 / 3 , leads to a Lagrangian that contains the 
terms: 


^ — (Xuus)ij '^R,i ‘S'e,1,4/3 + (Xeips)i '06,1,1/3 ^R,i ‘S'e,1,4/3 

+ {yQips)i Qi '06,1,1/3 • (5'3,2,1/6)^ + {Yids)ij dnj Tj • S' 3 ^ 2 , 1/6 + H.C.. (25) 


Naively, the 2-loop d = 7 contribution becomes more important than the 3-loop contribution, when the 
tree-level coupling of the necessary interaction TlriPs,2,7/677"^ is larger than Tu«sfQ,^s/(167r^) x (some 
logarithmic factor). 
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As above, we use the notation 1 , 4/3 = {Se^i, 4 / 3 )x{Te)fj and 1 , 1/3 = (V' 6 ,i,i/ 3 )x(T 6 )fj. 
Together, the Yq^ps, and Y^ds terms necessarily violate lepton number by two units. 

All generation indices in the couplings in Eq. (|25|) have been suppressed for simplicity. The 
contribution to the neutrino mass matrix can be roughly estimated as 


- 


Nc 


—— —(Yuus)33{Ye^s)a{YQp^s)3{YLds)30 + (« •H- /3) 
^'■LNV 


(26) 


where Alnv — 1 4/3 — 2 i/e — mass scale of the heavy states, which is 

typically taken to be TeV. Nf, is a colour factor. Here, we assumed that all the SM fermion 
masses are much smaller than Alnv. Putting all the Yukawa couplings in Eq. fl26]) equal to 
unity and Alnv = 1 TeV and Nc = 6 (for a colour sextet combination), one finds® 


{m^)ce ~ 1 X 10 ^ eV, ~ 2 x 10 ^ eV, {m^)rr ~ 0.3 eV. 


(27) 


This implies that the mass mechanism contribution to 0z//3/3 is guaranteed to be sub¬ 
dominant in this class of models. Also, Eq. fl271) shows that 3-loop models can potentially 
explain neutrino oscillation data only if all of the involved Yukawa couplings are set to be 
0{1). Thus, we expect such models to be quite constrained from upper limits on flavour 
violating decays of charged leptons. We will not discuss this class of models in more detail 
here, since their detailed phenomenology is outside the scope of this paper. 

The effective d = 9 operator resulting from the Lagrangian Eq. fl25l) can be written 
with the following linear combination of the basis operators of the short-range 

contributions to Offd/d decay as: 


C 


eff 


{Yuus)u{Yetps)e{YQp;s)l(YLds)el 




1 1 


— - —(05){|R}i? 


and the experimental bound Eq. 


16r " ^ 16 

constrains a combination of the coefficients to be: 


(28) 


{Yuus)ll{Ye'P;s)e{YQijs)l{YLds)el < 1 - 5'10 ^ 


y 1.0 [TeV^] ) ■ 


(29) 


The difference in the short-range bounds, Eq. (|2^ and Eq. ([29]), is due to the different values 
of nuclear matrix elements entering the transition operator. All other three-loop models will 
have constraints similar to the ones discussed here. They are listed in Table jVT] in the 
appendix. 
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Using (to;/) < 0.3 eV, we can formally write the constraint on the Yukawa couplings in the form of: 

(W„s)33(V^s)e(UQ,^s)3(Wds)e3 < 3 ■ 10^ , 

which is much worse than even the trivial constraint derived from perturbativity. 
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FIG. 11: An example of a four-loop neutrino mass model. To the left: Oi^/3/3 decay via the d = 9 
operator (ur UR)(d_R)(d_R)(FR Fr). To the right: four-loop d = 9 neutrino mass, see text. 


F. 4-loop models 


Finally, all operators = cr cr WRciRtrRciR, with exception of decomposition T-I-5-i 
(see Table H] in the appendix), will lead to four-loop neutrino mass models. The simplest 
possibility to construct a four-loop diagram for these operators is to use a SM charged 
current interaction. We estimate that this gives the dominant contribution to the neutrino 
mass. Here we show an example of the decompositions of the 0z//5/9 decay operator O- in 
Fig. [m which is based on decomposition T-I-3-ii (mr RR)((ij^)((i^)(eR cr). The four-loop 
neutrino mass diagram based on this decomposition is also shown on the right. Taking the 
limit ~ ~ ^ one can estimate the order of magnitude of 

this four-loop diagram, which is, 

~ 7Tr^2\A 


me^mepmuimujmdimdj 


(IGtt' 


^ 3 , 1 , 5/3 ‘^ 6 , 1 , 4/3 *^ 1 , 1,2 


i{^ees)afi^^uus)ij i^ipds')j- 


(30) 


The expression Eq. fl5U]) shows that this four-loop contribution would yield only ~ 

0(10“^°) eV for ~ uisq 14/3 ~ ’^ 51 , 1,2 ~ 1 TeV, even when choosing all SM fermion 

masses to be third generation. Since this is obviously many orders of magnitude below the 
values of neutrino masses required from oscillation experiments, models of this category 
by themselves cannot be considered realistic. Of course, neutrinos could be quasi-Dirac 
particles, explaining oscillation data by Dirac mass terms (using additionally introduced 
right-handed neutrinos), while Or/ 3/3 decay is dominated by the short-range diagrams such as 
the one shown in Fig. [HI However, constraints on Yukawa couplings will be similar to those 
derived in the previous subsections in Eq. flHll and Eq. fl28ll . with the exact value depending 
on the decomposition under consideration. All four-loop cases are listed in Table IVIII in the 
appendix. 


III. A CONCRETE 2-LOOP EXAMPLE 

In this section we will discuss one concrete genuine 2-loop neutrino mass model in some 
more detail. The example we choose is based on the decomposition T-I-4-ii-b of the Babu- 
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FIG. 12: Two loop diagram for neutrino masses generated by the Lagrangian in Eq. 


Leung operator On, which has not been discussed in the literature before. However, all Of(3(3 
decompositions that generate 2-loop neutrino masses behave quite similarly, in what concerns 
fits for neutrino oscillation data and constraints from lepton flavour violation searches. Thus, 
most of the discussion presented below can be applied qualitatively also to all other 2-loop 
decompositions. 

Any model of neutrino mass must not only generate the correct neutrino mass scale, 
but also be able to explain the observed neutrino mixing angles. For a recent update of 


all oscillation data, see, for example, 8^. In addition, since the neutrino mass matrix has 


a non-trivial flavour pattern, one also expects that low-energy models^® of neutrino mass 
are constrained by charged lepton flavour violation (LFV) searches. Here we will discuss 
only fi ey, since the experimental upper limit on this process provides usually the most 


stringent constraints in many models. We note that the authors of 83|] present a 2-loop 


model, which corresponds to the decomposition T-I-5-i and discuss also the constraints from 
other LFV searches, which we expect are very similar in our example. 

Below we will discuss two variations of the model based on T-I-4-ii-b. First (in Sec. IHI B|) . 
we introduce only one copy of the exotic fermion '06,2,-i/6 for simplicity. Next (in Sec. IIII C|) . 
we will allow to have three copies of these fermions, which allows to fit also quasi-degenerate 
neutrinos. 


A. General formulas for neutrino masses and p —> ey 

The Yukawa part of the Lagrangian describing the interactions between the exotic di¬ 
quark, S' 6 , 3 , 1 / 3 , the leptoquark, S' 3 pq/ 6 , and the coloured vector-like fermion, can be 

written as: 

(XQQs)ij{QiT ■ Q''j)SG^3,1/3 + (YL,ps)ak{La3^'ipk)Sls^^/^ (31) 

+ (kpds)fci('0fc<^iJ,i)*S'3,2,i/6 + {YliIs) ai 2,1/6 

10 By “low-energy” we mean TeV-scale, as in contrast to “high-scale” seesaw models. 


21 















Here, i,j are generation indices for quarks, we use Greek indices for lepton generations and 
k runs over the number of copies of ' 06 , 2 ,-i/ 6 - This Lagrangian generates a 2-loop diagram 

Following the general 


63 


which corresponds to PTBM-4 according to classihcation by 
formulas from 63j], the neutrino mass matrix can be expressed as: 

^ [{yQQs)ij{^Lips)ak{yiljds) kiO^Lds)j3j T 0^QQs)ij0^Lips)f3k0^ipds)ki0^Lds)aj] 


X , rud ,, ^56,3,1/3, ^ 


(32) 


where Nc is a colour factor, with Nc = Q for this model. Summation over all flavour indices 
i,j,k is implied. ^ ^/g, rurf., 3 ^^ 3 , is a loop integral dehned as: 


X / (fq / (fk 


F(jn^i^ , '^832,1/6 ’ ’ ^'S'6,3,1/3 ’ ) 

1 


mdiTUd. 


vr^ 


(33) 


{(f - ml )(g 2 - )(fc 2 - ml){k^ - ){{q + ky - ) 


Due to the strong hierarchy in down-type quark masses, the integral in Eq. fl3^ is completely 
dominated by the contributions from bottom quarks, unless the couplings Fg-gg, and 
'^Lds follow an equally strong inverse hierarchy. We have thus taken into account only the 
contributions from bottom quark exchange in our numerical evaluation. Since it is convenient 
to rewrite Eq. fl3^ in terms of dimensionless parameters, we dehne 2 , r and tk as 


rriq 

*3,2,1/6 

.2 = - 5 -. 

mt 


rriq 
*6,3,1/3 

r = - 

mi 


and 


tk — 


_ 


Rescaling the loop momenta, the integral can then be written as: 


mt 


(34) 


1 ( 4 , 2 ;, l,r) = 


1 




/ d‘^k 


(d^ — 4)(<?^ — z){k?‘ — l)(fc2 — r)((g -I- ky — 1) 


(35) 


This integral has been analytically calculated several times in literature. We follow the 
procedure outlined in j^, based on the calculations of js^. We will ht the neutrino mass 
calculated with Eq. fl32l) to neutrino oscillation data. The discussion depends on the number 
of copies of the fermion mediator 46 , 2 ,-i/ 6 ; as mentioned above we will discuss two different 
scenarios in the following subsections. 

The rate of the LEV process p —)■ ey has also been calculated several times in literature. 
We adapt the general formulas shown in 85j for our particular case. The amplitude for 
/i —)■ ey decay is given by 


A4(p -)■ ey) = eaRe*^qf}u{pe)i(j'^^u{p^), (36) 

where e is the electric charge, is the photon polarization vector, q^ is the momentum 
of photon, and cr"^ = (i/2)[y",y^]. There are two contributions to the coefficient in 
the model we are discussing; one is the one-loop diagram with the diquark and the exotic 


22 















fermion, the other is that with a bottom quark and the leptoquark. The total an is given 
by 


. rn. 


o-R = I 






\ ^ ^2i^k) I ^ N, 2 T2 (^s) -^l(^s) 

L'lpSJke^^ r fr MSj/xSf-'M5/3e- 


mt 


56,3,1/3 


mt 


53,2,1/6 


. (37) 


where Xk = 2 — and xs = 2 '"^ —• The functions Fi{x) and F 2 {x) are defined as 

®6,3,l/3 ‘®3,2,l/6 


Fi{x) 

F 2 {x) 


— 5x — 2 X Inx 


+ 


12(x-l)3 2(x-l)4’ 

2x^ + 5x — 1 In X 


12(x-l)3 2(x-l)4’ 


(38) 

(39) 


which are presented in Eqs. (40) and (41) in 8^. The branching ratio for the p ^ ey 
process, neglecting the electron mass, can then be expressed with the coefficient an as 


Br(/i —)■ ey) ~ 


487r^Q! IctrI^ 

Glml ’ 


(40) 


where a is the fine-structure constant. 


B. One generation of V’ 6 , 2 , 1/6 


The analysis presented in this section uses very similar methods to the one in ref. 16| . 
where double beta decay and LFV is discussed in a 1-loop neutrino mass model containing 
colour octets. We will first consider a variant of the model, in which there is only one copy 
of the fermion mediator '06,2,1/6- The expression for the neutrino mass matrix in this case is 
given by suppressing the index for '0 in Eq. fl32l) . which gives 


where 


{'m„)ai3 = [{yLips)a{yLds)l33+ {yLips)p{yLds)a3]J^, 
Njn^ 


F = 


:{yQQs)33{ytpds)3l I;''")- 


(41) 


(42) 


(Ibvr^)^ 

Since det(mi,) = 0 in this case, this version of the model can fit only to the hierarchical 
neutrino mass spectra (both of the normal and the inverse type), but not to the degenerate 
spectrum.The eigenvalues of Eq. flTT]) can be easily found to be: 


m 


*31(3) 


= 0 , m 


7^2,3(1,2) 


a{yLds)c:i 1 = a/ l(U,/,5)ct ^ l(lLrfs)u 3 | 


T (43) 


However, we remind that this is true, only when contributions to neutrino masses from the first and the 
second generation quarks are negligible. 
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Meir(TeV) 


FIG. 13 : The prefactor T^ defined in Eq. (I 4 ^ . for {Yqqs) 33 = {Y^ds)3 = 1 in units of keV as a 
function of in TeV. Here, M^g = rn^ ~ ms^ ^ i/e — 31/3- 

for normal hierarchy (inverted hierarchy). In Fig. [131 we give typical values for the common 
factor T, which are calculated with the assumption of a nearly degenerate spectrum of heavy 
particles with the mass scale M^g = ~ "^ 53 , 2 , 1/6 - ^Se, 3 , 1/3 and (Yqqs )33 = (Xipds )3 = 1- 

From Eq. (H3|l and Fig. [131 one can estimate the constraints from neutrino masses on the size 
of the Yukawa couplings. In order to reproduce the neutrino mass suggested by atmospheric 
neutrino oscillation ~ 0.05 eV), keeping the common mass scale M^g at 1 TeV, the 
Yukawa couplings Y£:^g and Y^^g must be set typically to 0(10“^). 


Although the eigenvectors of Eq. fjTTl) can be calculated analytically, numerical exercises 
might be more helpful to grasp phenomenological aspects of the model. In the following, 
we will generate random sets of Yukawa couplings (Yr 7 />.'?) »■ an d {YLds)i 33 under the condition 
that they reproduce the latest neutrino oscillation data within 3 a C.L. We will only 
show plots with the Yukawa couplings that £t the normal hierarchical neutrino spectrum, 
because plots for the inverse hierarchical case look qualitatively similar. 

Let us start the discussion with double beta decay. The half-life of 0^(3(3 induced by the 
Majorana mass of neutrino is proportional to the inverse-square of the effective neutrino 


mass: 


yiOupfi 
■ 1/2 


OC 


[m„ 


-2 


(44) 




For the normal hierarchy case, the effective mass is roughly given as ~ 5 ^ 2 " 

3 X 10“^ eV, which results in half-lives of the order of ~ 10^® ys. For the inverse 

hierarchy case, one Ends {m^)ee i/Am^ ~ 5 X 10 ^ eV, which leads to Ty 2 ^ ~ 10^'^ 
ys. The current experimental limits to the half-life of ^^®Xe and ^®Ge are of the order of 

rynOL'/S/? 

1/2 

rjiOu^jS 
1/2 

lifes. 


(1-2) X 10^^ ys 0, 11,18^, while the next round of experiments could reach eventually 
10^^ ys. Therefore, only the inverse hierarchical case can result in measurable half- 
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Meff(GeV) 


FIG. 14: Calculated half-lives for 0vl3l3 decay of ^^®Xe, considering only the short range contri¬ 
bution to the decay rate. The different colours correspond to (from left to right) = 1, 5, 10 
and 50. If the third and first generation couplings are of the same order, Oz//3/3 decay will have an 
immeasurably large half-life in the variant of the model with only one copy of ' 06 , 2 , 1 / 6 - 


The short range-contribution due to the d = 9 operator (cf. Eq. (| 2 T|) ) is proportional to the 
following combinations of the parameters: 


1/2 


OC 


(XqQs) 11 e (Ei/ds) 1 iyLds)el 


n -2 




(45) 


i.e., while the neutrino mass matrix is dominated by Yukawa couplings of the third quark 
generation, double beta decay is sensitive only to the Yukawa couplings that couple to the 
hrst generation quarks. To discuss the relation between these two contributions to Oufdfd, 
we introduce a scaling factor 


h31 


(XQQs)ll(Xiljds)l(XLds)el 

( Eqqs) 33 ( ) 3 ( ) e3 



(46) 


i.e., 7/31 = 1 corresponds to quark flavour universality in the Yukawa couplings. In Fig. [HI 
we calculate half-lives induced from the short-range contribution with randomly generated 
Yukawa couplings, assuming different values of r /31 G {1,5,10,50}. Taking 7/31 = 1, we hnd 
quite long half-lives, too large to be measured in realistic experiments. On the other hand, 
with 7/31 = 10, we hnd a lower limit on which is approximately > 400 GeV. This is 
still not competitive with leptoquark searches at the LHC, which places constraints on the 
masses of leptoquarks at 77753 ~ (^00 ~ 1000) GeV (depending on generation) already 


in the hrst run [87H 89| Thus, it is reasonable to conclude that, as for the mass mechanism 


contribution, the short-range contribution to the half-life is also expected to be too long to 
be measured in the near future in this variant of the model, unless 7/31 is very large (i.e., for 
highly inverse hierarchical Yukawa couplings in terms of the quark generations). 
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FIG. 15: Branching ratio Br(/r —>■ 67 ) as a function of Meg in GeV. Red (blue) points have been 
calculated with (Fqqs )33 = ( 1 ^m ,‘?)3 = 10“^ (10“^). The horizontal line is the experimental upper 
limit from the MEG experiment 9C|. 


Finally, we discuss briefly the LFV process /i ey. While the neutrino mass matrix is 
proportional to the combination of the Yukawa couplings (YLp5)a(kLds)/33(^QQs)33(^V'rf5)3, 
the branching ratio of Br(/i -)■ ey) depends only on \{YL^s) 2 {yl^s)i\^ |(YLd5)23(^2d5)i3p- 

In Fig. [m we show Br(/i ey) for two different choices of the set of {Yqqs)33 and {Y^cis)3, 
as a function of M^s, assuming again that the mass spectra of heavy particles are nearly 
degenerate for simplicity. With the choice {Yqqs) 33 = {Y^cis)3 = 10“^ the LFV process can 
place a bound on of roughly Mes > TeV. However, the bound depends strongly on the 
exact choice of the remaining Yukawa couplings and Yids- On fhe other hand, the LFV 
process can exclude only few parameter points in the case of (Yqqs) 33 = {Y^ds)3 = 10“\ 
and no useful limit on Mefj can be derived. 

As we have seen in this subsection, this variant of the model can reproduce oscillation 
data without running into conflict with LFV searches. However, it is interesting to note that 
the {Yqqs)33 interaction with the size required for reproducing neutrino masses will result 
in sizeable decay rates of the diquark into third generation quarks (both tops and bottoms), 
which should be testable at the LHC. 


C. Three generations of V’6,2,1/6 

Next, we examine the model with more than one copy of the fermion mediator "06,2,-1/6) 
which can fit not only hierarchical neutrino mass spectra, but also a quasi-degenerate spec¬ 
trum. Here, we introduce three copies of "06,2,-1/6) motivated by the observed generations of 
SM fermions. 

To simplify the following discussion, we adopt the following ansatz in the flavour structure 
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of the Yukawa couplings: 


(XLds)a 3 {y^xpds)k 3 — y(XLips)ak- 


( 47 ) 


With this ansatz, all the flavour structure relevant to phenomenology can be represented 
with only one vector (apart from a possible normalization factor y). The neutrino mass 
matrix can then be cast into the form: 




where the A is dehned as 


and I is given as 


AqA: = iXL'tl}S)ak — —iXLds)a3(XTpds)k3-, 


h — - 7 T 7 r^^y{^QQs) 33 l{z, 1 , 4 y) 


( 48 ) 


( 49 ) 


(IGtt' 


( 50 ) 


Comparing Eq. fl 48 |) with the neutrino mass and mixing matrix, we can hnd the direct 
relation between A and the measured neutrino data. Following the procedure originally 
developed by Casas and Ibarra for seesaw type-I | 9 l|, we parametrize A as 




( 51 ) 


Here, rhy is the matrix of eigenvalues of rui,, which is diagonalized with the neutrino mixing 
matrix Up via 

(Up ')ia {j^v)ay (,Up')yj = tflp diag ^TTlp^ 'kTlp2 TTlpJ^ , (^2) 

for which we use the following standard parametrization 


( 53 ) 


Cij = cos 6ij, Sij = sin6*jj with the mixing angles 5 is the Dirac phase and ai, a2 are 
Majorana phases. Finally, i? is a complex orthogonal matrix which satishes the condition 
R = 1 . We use the following parametrization for the R matrix in terms of three complex 
angles 6*i, 6*2, and 63 as 



( 

C12C13 

S12C13 

5136*"^ ^ 


/ giai 0 

0^ 

Up = 


-S12C23 - Ci25235i3e"*^ 

C12C23 - 5125235136"*^ 

523C13 


0 6*'^^ 

0 


V 

S12S23 - Ci2C23Sl3e“*'^ 

-612523 - 5126235136"*'^ 

C23C13 j 


0 

0 

V 


R = 


( C2C3 —C1S3 — S1S2C3 S1S3 — C1S2C3 \ 
C2S3 C1C3 — 515253 —51C3 — Ci5253 

y 52 5 iC2 C 1C2 j 


( 54 ) 


12 


This ansatz can be justified by introducing a flavour symmetry with flavour-charged scalar (flavon) fields. 
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After fitting the neutrino oscillation data with the parametrization shown above, there 
remain y, (Yqqs )33 and the masses as free parameters. For the 

calculation of the short-range contribution to the 0z//5/9 decay, we also have the parameter 
r] 3 i. For simplicity, we set y = 1 and assume again a nearly degenerate spectrum for heavy 
particles, which is parameterized with Meg. We can then calculate half-lives for both 

neutrino mass mechanism and the short-range contribution, as a function of M^s and 
ri 3 i. In Fig. [161 we hxed the oscillation parameters and at their best-£t 

values, while s^g = 1/2 and 5^2 = 1/3 and set 6 as well as the Majorana phases to zero, just 
to sketch out some phenomenological aspects of this example. Each panel shows the half-life 
of the short-range contribution for decay of ^^®Xe as a function of (top panel), Meg 
(middle panel) and (bottom panel). In each panel, we examine several choices for the 
remaining parameters, which are explained in the hgure caption. The corresponding half- 
lives induced from the mass mechanism are also indicated. As shown in Fig. [161 half-lives 
can vary over many orders of magnitude with the choice of parameters. The amplitudes 
induced from the mass mechanism becomes the same order as that from the short-range, 
when ri 3 i ~ 2.7 (6.5) for Meg = 0.5 TeV (1 TeV). As in the case with only one generation of 
t/ 6 , 2 , 1 / 6 , the mass mechanism dominates the 0z//3/3, if the ratio rjsi is taken to be unity and 
the heavy mass scale Meg is given at the typical LHC search sensitivities. However, since the 
three-generation case can £t the quasi-degenerate neutrino spectrum, decay half-lives 

can be much shorter than in the one generation case and can saturate the experimental 
bound. 

We now turn to Br(p —)■ ey). Again, as in the one generation case, the neutrino mass 
matrix depends on Yukawa couplings, but is not directly related to Br(/i —)■ ey). Therefore, 
we have always the freedom to adjust {Y^ds)k3 and (Yqqs)33 so as to fit the neutrino masses. 
The other Yukawa couplings are then hxed by the neutrino data (and the choice of Meg), 
and we can use them to calculate Br(/i ey). Fig. [T7I shows some examples with a value of 
Meg = 1 TeV. The plots show that constraints from Br(p —)> ey) can be easily fulhlled. For 
this choice of Meg, only if both {Yqqs)33 and {Y^ds)i 3 are set to order 0(10“^) or lower, the 
predicted Br(p —>■ ey) can saturate the experimental bound. 


IV. CONCLUSIONS AND DISCUSSION 

We have discussed the relation between the d = 9 short-range contributions to the OujSfd 
decay amplitude with neutrino mass models. All contributions to 0i^/3/3 decay violate lepton 
number and, therefore, generate also Majorana neutrino masses. We have classihed all possi¬ 
ble (scalar-mediated) short-range contributions to the decay rate according to the loop level, 
at which the corresponding models will generate Majorana neutrino masses. Possibilities 
range from tree-level to 4-loop neutrino masses. For each case we have discussed one exam¬ 
ple briehy and given estimates of the typical constraints imposed by both the short-range 
contribution and the mass mechanism. Generally, one expects that for models with tree- 
or 1-loop neutrino masses, the short-range Offd/d decay amplitude will be sub-dominant to 





0.20 0.30 0.50 0.70 1.00 1.50 2.00 3.00 


Mrff(TeV) 




Merr(TeV) 



FIG. 16: Calculated half-lives for 0uf5f5 decay of ^^®Xe, considering only the short-range contri¬ 
bution to the decay rate. The various plots show from top to bottom: T 1/2 versus 7731 , and 
rriui, for a fixed set of neutrino oscillation parameters and different choices of 7731 , Mes and rriy^ as 
follows: In the top plots to the left (right) nii,^ = 0.3 eV (0.05 eV), different lines show different 
choices of Mgfj; from left to right: Meg = 0.2, 0.5, 1, 2 and 5 TeV. In the middle panel, to the 
left (right) = 0.3 eV (0.05 eV), different lines show different choices of 7731 ; from left to right: 
7731 = 1, 2, 5, 10 and 20. In the lower panel, to the left (right): Meg = 0.5 TeV (1 TeV), different 
lines are for different choices for 7731 ; from top to bottom: 7731 = 2, 3, 5, 8 and 10. For comparison 
we also show the half-lives for the neutrino mass mechanism as horizontal lines in the top and 
middle panel and as cyan lines in the lower panel. Oscillation parameters are sfg, and ^rn^i 

at their best-fit values, while S 23 = 1/2 and 5^2 = 1/3 for the case of normal hierarchy. 
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FIG. 17: Br(/i e'y) versus Y^ds = iYipds)i 3 (left) and Yqqs = (1qq5)33 (right), for a fixed choice 
of Mgff = 1 TeV, neutrino oscillation parameters as in Fig. [16] and mi,^ =0.05 eV (red lines) and 
rriui — O-H eV (blue lines). Full, dot-dashed and dashed lines are for Y^ds (right) and Yqqs (left) 
equal to 1, 10“^ and 10“^ respectively. Br(^ —>■ ey) can saturate the experimental bound only 
for small values of these couplings, since smaller values of Y^dS arid Yqqs require larger values for 
Yi^ps and hr order to fit neutrino data. 

the mass mechanism. For 2-loop models short-range Oir/3/3 decay amplitude and mass mech¬ 
anism can be comparable, while for 3-loop and 4-loop models the short-range part of the 
amplitude will dominate. 

We have also discussed one particular example of a 2-loop model in more detail. Here, 
we have shown different parts of parameter space where mass mechanism or short-range 
amplitude dominant can each be dominant. In the study, we have taken recent neutrino 
oscillation data and constraints from LFV experiments into consideration. 

In the appendix we give the full list of decompositions, classihed according to our scheme, 
in tabular form. 
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V. APPENDIX 

Here, we give tables in which all possible scalar short-range decompositions are classihed 
according to the loop level at which they will generate neutrino masses. The decompositions 
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which generate neutrino masses at tree, 1-loop, 2-loop, 3-loop and 4-loop level are listed in 
Tables ID to I VIII The identification number given to each decomposition is defined in 3l|, . 

The notation T-I and T-II refers to the two possible topologies of the decompositions of 
d = 9 Oufdld effective operators, and the BL number is the Babu-Leung classification of the 


effective neutrino mass operator, given in 36|]. The columns “Add. Int.” specify additional 
interactions, with respect to those appearing in the decomposition. While they do not 
appear directly in the Oz//3/3 diagram, these additional interaction can not be forbidden by 
any symmetry, without forbidding the corresponding 0i^/3/3 decay decomposition at the same 
time. Once present, they generate a neutrino mass diagrams at the quoted loop level. The 
columns “Diagram” specify the topology of the neutrino mass diagram, for example “type- 
I” for seesaw type-I and so forth. The identification numbers for the 1-loop and 2-loop 
neutrino mass diagrams are taken from the general topology classification given in 78|] and 


63| |. respectively. 


Here, we briefly comment on “associated operators”. As discussed in Ref. 78|], some of the 
decompositions generate not only the original effective operator but also necessarily generate 
other operators, when all possible contractions are carried out. We call this associated 
operators. For example, the decomposition T-I-2-iii-a of the BL ^19 operator consists of 
the following fundamental interactions. 


/lT-I-2-iii-a —yidsiLdu) ■ 5'3,2-1/6 + ' >51 2 _i/6 

+ 1 , 1/3 (yXdn) + Yues >5g, 1 , 1/3 + H.c., 


(55) 


where A is the Gell-Mann matrices. The first two interactions, together with the Majorana 
mass of the fermion '08,i,o result in the BL ^11 operator {LdR){Q){Q){LdR). In the same 
way, the last two interactions lead to the d = 9 lepton number violating effective operator 
{uReR) {dR) (dR) (urCr) , which is in Eq. ([5]). ^1 the decompositions accompanied by 


associated operators were listed in tables of Ref. [78| . We take into account the associated 


effective operators in our classification scheme. In short, if the associated operator generates 
neutrino masses at a lower loop level than the original one, we classify the decomposition 
with the loop level of the associated operator. The Lagrangian for a concrete example is 
given in Eq. fl55ll . Here, although the original effective operator BL ^19 gives neutrino 
masses only at the 3-loop level, the decomposition T-I-2-iii-b of BL ^19 also produces BL 
and it generates the 2-loop neutrino mass diagram with the help of the SM Yukawa 
interactions. Therefore, we list the decomposition T-I-2-iii-b of BL ^19 as a 2-loop neutrino 
mass model in Tab. lYl More examples are given in the tables. 
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FIG. 18; Dimension 7 {d = 7) neutrino mass diagrams generated by the decompositions listed in 
Table El 
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T-I # Op. 

BL# 

S 


S' 

Diagram 

Add. Int. 

1-i {ud){e){e){ud) 

11, 12, 14 


(l,l)o 

(1,2)i/2 

type I 

LfiuoH^ 

1-i {ud){e){e){ud) 

11, 12, 14 

(1,2)+i/2 

(l,3)o 

(1,2)i/2 

type III 


1-ii-a {ud){u){d){ee) 

11, 14 

(1>2)-|-1/2 (3, 3)+2/3 

(l,3)+i 

type II 


1-ii-a {ud){u){d){ee) 

11, 14 

(8, 2)-hi/2 (3,3)+2/3 

(l,3)+i 

type II 


1-ii-a {ud){u){d){ee) 

12, 14 

(1)2),i/ 2 (3, 2),7/g 

(l,3)+i 

type II 


1-ii-a {ud){u){d){ee) 

12, 14 

(8, 2)+i/2 (3, 2)+7/6 

(l,3)+i 

type II 


1-ii-b {ud){d){u){ee) 

12, 14 

(1)2),7/2 (3, 3),i/3 

(l,3)+i 

type II 


1-ii-b {ud){d){u){ee) 

12, 14 

(8, 2)+i/2 (3, 3)+i/3 

(l,3)+i 

type II 


1-ii-b {ud){d){u){ee) 

11, 14 

(1)2)-|-1/2 (3, 2),5/3 

(l,3)+i 

type II 


1-ii-b {ud){d){u){ee) 

11, 14 

(8, 2)+i/2 (3, 2),5/3 

(l,3)+i 

type II 


2-i-b {ud){e){d){ue) 11, 19, 14, 20 

(1,2)+i/2 

(l,l)o 

(3,1)+i/3 

type I 


2-i-b {ud){e){d){ue) 

11, 14 

(1,2)+i/2 

(l,3)o 

(3, 3)+i/3 

type III 


2-ii-b {ud){e){u){de) 

11, 14 

2),i/2 

(l,l)o 

(3, 2)+i/6 

type I 


2-ii-b {ud){e){u){de) 

11, 14 

(1,2)+i/2 

(l,3)o 

(3, 2)+i/6 

type III 


2-iii-a {de){u){d){ue) 

11, 19 

(3,2)-i/6 

(l,l)o 

(3,1)+i/3 

type I 

LfiuoH^ 

2-iii-a {de){u){d){ue) 

11 

(3,2)-i/6 

(l,3)o 

(3, 3)+i/3 

type III 


3-ii {uu)[d){d){ee) 

11 

(6, 3)+i/3 (3,3)+2/3 

(l,3)+i 

type II 


3-ii {uu){d){d){ee) 

12 

(6, l)+4/3 (3, 2)+7/6 

(l,3)+i 

type II 


3-iii {dd){u){u){ee) 

12 

(6,3)i/3 

(3, 3)+i/3 

(l,3)+i 

type II 


3-iii {dd){u){u){ee) 

11 

(6, l)+2/3 (3, 2)+5/6 

(l,3)+i 

type II 


4-i {de){u){u){de) 

11 

(3, 2)i/6 

(l,l)o 

(3, 2)+i/6 

type I 

LfiuoH^ 

4-i {de){u){u){de) 

11 

(3,2)-i/6 

(l,3)o 

(3, 2)+i/6 

type III 

Lijji3oH^ 

5-i {ue){d){d){ue) 

11, 19, - 

(3,1)-i/3 

(l,l)o 

(3, l)+l/3 

type I 

LfinoW 

5-i {ue){d){d){ue) 

11 

(3,3)i/3 

(l,3)o 

(3, 3)+i/3 

type III 



d, 

O 

1 —1 

1 —1 

BL# 

S S' S" 

Diagram 

Add. Int. 

1 

{ud){ud){ee) 

11, 12, 14 (1,2)+i/2 (1,2)+i/2 (1,3)i 

type II 

S 13 - 1 HH 

1 

{ud){ud){ee) 

11, 12, 14 (8,2)+i/2 (8,2)+i/2 (1,3)i 

type II 

S 13 - 1 HH 

3 

{uu){dd){ee) 

11 

(6,3)+i/3 (6,1)+2/3 (1,3)-i 

type II 

S 13 - 1 HH 

3 

{uu){dd){ee) 

12 

(6,1)+4/3 (6,3)i/3 (1,3)i 

type II 

S 13 - 1 HH 


TABLE I: List of the decompositions that generate neutrino masses at tree level. The ID-numbers 

with “T” are assigned as in Ref. iSil . and the decomposition is specified in the “Op. ” column. We 

also give the ID-numbers of lepton-number-violating effective operators, which are classified as in 

Babu and Leung iSa l. in “BL#”. The SM charges of fields appearing in the decomposition are also 

given. “Diagram” indicates the type of resulting tree-level neutrino mass diagrams: “type I” for 

type I seesaw mechanism, and so on. In the column “Add. Int. ”, we give the additional interaction 

that is missing in the decomposition but is necessary to generate the neutrino mass diagram. For 

the decompositions in this table, unless some severe fine-tuning of parameters is done, the mass 

mechanism of double beta decay will dominate over the short-range contributions. 
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H 

1 

1—1 

O 

T 

BL# 

5 

if S' 

Diagram 

Add. Int. 

1-i {ud){e){e){ud) 

11, 12, 14 

(8>2),i/2 

( 8 ,l)o ( 8 , 2 )i /2 

Ti/-3 


1-i {ud){e){e){ud) 

11, 12, 14 

(8) ‘^)+l/2 

(8,3)o (8,2)3/2 

Ti/-3 


2-i-a {ud){d){e){ue) 

11, 14 

(3, 2),5/6 (3, l),i/3 

Ti/-l-iii 


2-i-a {ud){d){e){ue) 

11, 14 

(1)2),i/2 (3,2),5/6 (3,3)+i/3 

Tz^-l-iii 

dR'^32^H] 

2-i-a {ud){d){e){ue) 

11, 14 

(8, 2)+i/ 2 (3,2),5/6 (3, l)+i/3 

Ti/-l-iii 


2-i-a {ud){d){e){ue) 

11, 14 

(8, 2),i/2 (3,2),5/6 (3,3)|i/3 

Tu-l-iii 


2-i-b {ud){e){d){ue) 11, 14, 19, 20 

(8, 2),i/2 

(8,l)o (3,1 )+i/3 

Ti/-3 


2-i-b {ud){e){d){ue) 

11, 14 

(8, 2)+i/2 

(8,3)o (3,3 )|i/3 

Ti/-3 

S82lS82.^^t^^t 

2-ii-a {ud){u){e){de) 

11, 14 

(1)2),i/2 (3,3)+2/3 (3, 2),i/6 

Tz^-l-iii 


2-ii-a {ud){u){e){de) 

11, 14 

(8, 2)+i/2 (3,3)+2/3 (3, 2)+i/6 

Tz^-l-iii 


2-ii-b {ud){e){u){de) 

11, 14 

(8, 2)+i/2 

(8,l)o (3,2 )+i/6 

Ti/-3 

^82l^^82fdd^ 

2-ii-b {ud){e){u){de) 

11, 14 

(8, 2)+i/2 

(8,3)o (3, 2 ),i/6 

Ti/-3 


2-iii-a {de){u){d){ue) 

11 

(3,2)_i/6 

(8,l)o (3,1 )+i/3 

Tu-l-ii 


2-iii-a {de){u){d){ue) 

11 

(3,2)_i/6 

(8,3)o (3,3)_|_3/3 

Tu-l-ii 


2-iii-a {de){u){d){ue) 

14 

(3,2)_i/6 

(1,2)+3/2 (3,1)+i/3 

Tu-l-ii 


2-iii-a {de){u){d){ue) 

14 

(3,2)_i/6 

(8, 2)+i/2 (3, 1)+i/3 

Tu-l-ii 


2-iii-a {de){u){d){ue) 

14 

(3,2)_i/6 

(1>2)+i/2 (3,3)+i/3 

Tu-l-ii 

4-iSs3iH' 

2-iii-a {de){u){d){ue) 

14 

(3,2)_i/6 

(8, 2)+i/2 (3,3)+i/3 

Tu-l-ii 

4-A33iff' 

2-iii-b {de){d){u){ue) 

11 

(3,2)_i/6 

(3,2)_,_3/g (3, l)_,_i/3 

Tu-l-ii 

4 .S 3 . 1 ff' 

2-iii-b {de){d){u){ue) 

11 

(3,2)_i/6 

(6, 2)+i/6 (3, 1)+i/3 

Tu-l-ii 

4-iff3i|ff' 

2-iii-b {de){d){u){ue) 

11 

(3,2)_i/6 

(3, 2)+i/6 (3,3)+i/3 

Tu-l-ii 

4-4s3|ff' 

2-iii-b {de){d){u){ue) 

11 

(3,2)_i/6 

(6,2)_,_3/g (3,3 )_,_i/3 

Tu-l-ii 

4-433|ff' 

2-iii-b {de){d){u){ue) 

14 

(3,2)_i/6 

(3,1)_i/3 (3,1)+i/3 

Tu-l-ii 

4-43.iff' 

2-iii-b {de){d){u){ue) 

14 

(3,2)_i/6 

(6,1)_i/3 (3,1)+i/3 

Tu-l-ii 

4-43.iff' 

2-iii-b {de){d){u){ue) 

14 

(3,2)_i/6 

(3, 3)_i/3 (3,3)_|_i/3 

Tu-l-ii 

4-433iff' 

2-iii-b {de){d){u){ue) 

14 

(3,2)_i/6 

(6, 3)_i/3 (3,3)+i/3 

Tv-l-n 

4-433iff' 

3-i {uu){e){e){dd) 

11 

(6, 3)_,_i/ 3 (6, 2)_3/6 (6, l)_2/3 

Ti/-3 

S33iSL|ifH' 

3-i {uu){e){e){dd) 

12 

(6, l)+4/3 (6, 2)+5/6 (6, 3 )+i/3 

Ti/-3 

s..|4iff*ff’ 

4-ii-a {uu){d){e){de) 

11 

(6, 3)+i/3 (3,3)+2/3 (3, 2),i/6 

Tz^-l-iii 

Q''l'33^dd^ 

5-ii-b {ue){e){u){dd) 

11 

(3, l)-l /3 

(3, 2)_5/g (6, l)_2/3 

Ti/-l-iii 

4?332-|ff' 

5-ii-b {ue){e){u){dd) 

11 

(3,3)_i/3 

(3, 2)_5/g (6, l)_2/3 

Tz/-l-iii 



TABLE II: Decompositions that generate neutrino masses at 1-loop. The namin g c onvention of 
1-loop neutrino mass diagram, which is used in “Diagram” column, follows Ref. 178 ] and is also 


shown in fig. 0- For the decompositions in this table, unless some severe fine-tuning of parameters 

is done, the mass mechanism of double beta deeay will dominate over the short-range contributions. 
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T-II # Op. BL # 5 S' S" Diagram Add. Int. 

2 (nd)(ixe)((ie) 11, 14 (1,2 )_h/2 (3, l)_i/3 (3,2)_i/6 Tu-l-ii i i 

2 (nd)(ne)(de) 11, 14 (1,2 )+i/2 (3,3)_i/3 (3,2)_i/e Tu-l-ii ii 

2 (nd)(ne)(de) 11, 14 (8,2)+i/2 (3, l)-i/3 (3,2)_i/6 Tz/-l-ii i i 

2 {ud){ue){de) 11,14 (8,2)+i/2 (3,3)_i/3 (3,2)_i/6 Tz/-l-ii 5 't2_ 15 '^3_ i ^ 

6 3 


TABLE III: Decompositions (T-II) that generate neutrino mass at 1-loop, which are continued from 
Tab.m For the decompositions in this table, unless some severe fine-tuning of parameters is done, 
the mass mechanism of double beta decay will dominate over the short-range contributions. 


T-I # Op. 

BL# 

S 

# 

S' Diagram 

2-iii-a {de){u){d){ue) 

19 

(3,2)-i/6 

(8,l)o 

(3,1)+i/3 PTBM-1 

4-i {de){u){u){de) 

11 

(3,2)-i/6 

(8,l)o 

(3,2)+i/6 PTBM-1 

4-i {de){u){u){de) 

11 

(3,2)-i/6 

(8,3)o 

(3,2)+i/6 PTBM-1 

4-ii-b {uu){e){d){de) 

11 

(6) 3)+i/3 (6, 2)3/g 

(3,2)+i/6 PTBM-4 

5-i {ue){d){d){ue) 

11, 19 (3,1)i/3 

(8,l)o 

(3,1)+i/3 PTBM-1 

5-i {ue){d){d){ue) 

11 

(3,3)i/3 

(8,3)o 

(3,3)+i/3 PTBM-1 

5-ii-a {ue){u){e){dd) 

11 

(3) l)-l/3 

(6,2 )-i/6 

(6, l)2/3 PTBM-4 

5-ii-a {ue)[u){e){dd) 

11 

(3,3)i/3 

(6,2 )-i/6 

(6, l)2/3 PTBM-4 


d 

O 

1—1 

1—1 

BL# 

S 

S' 

S" Diagram 

4 {uu){de){de) 

11 

(6 , 3)|i/ 3 (3,2)i/g 

(3,2)i/6 CLBZ-1 

5 {ue){ue){dd) 

11 

(3,1)-i/3 

(3,1)-i/3 

(6) l)-i-2/3 CLBZ-1 

5 {ue){ue){dd) 

11 

(3,3)i/3 

(3,3)i/3 

(6 , 1)-i-2/3 CLBZ-1 


TABLE IV: Decompositions that generate the d = 5 neutrino mass operator LLHH at 2-Loop. We 
follow the naming convention used in Although the effective operator of BL #19 can generate 


neutrino mass only at the 3-loop level ]3n . \3't] . the decompositions of BL #19 listed in this table 


generate not only the BL #19 but also the “associated” BL #11 operator and thus are classified as 
2-loop neutrino mass models. For the decompositions in this table, the mass mechanism of double 
beta decay and the short-range contributions can he comparable. 
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[20] N. Steinbrink et al.. New J.Phys. 15, 113020 (2013), arXiv: 1308.0532. 

[21] J. Lesgourgues and S. Pastor, Phys.Rept. 429, 307 (2006), arXiv:astro-ph/0603494. 

[22] S. Hannestad, Prog.Part.Nucl.Phys. 65, 185 (2010), arXiv: 1007.0658. 

[23] Y. Y. Wong, Ann.Rev.Nucl.Part.Sci. 61, 69 (2011), arXiv:1111.1436. 
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T-I # Op. 

BL# 

S S' 

Diagram 

Add. Int. 

2-i-a {ud){d){e){ue) 

19, 20 (1,2)+i/2 (3,2)+5/6 (3,1)+i/3 

(b) 


2-i-a {ud){d){e){ue) 

19, 20 (8,2)+i/ 2 (3,2)+5/6 (3, l)+i/3 

(b) 


2-ii-a {ud){u){e){de) 

19, 20 

(3, 2)+7/6 (3, 2)+i/6 

(a) 


2-ii-a {ud){u){e){de) 

19, 20 

(3, 2)+i/2 (3, 2),7/g (3, 2),3/g 

(a) 


2-ii-b {ud){e){u){de) 

19, 20 

+ 0--,‘^)-l/2 (3, 2),l/g 

(c) 

eRipi2-lH^ 

2 

2-iii-a {de){u){d){ue) 

19 

(3,2)_i/g (8,l)o (3,1 )+i/3 

(d) 


2-iii-a {de){u){d){ue) 

20 

(3,2)_i/g (1,2 )+i/2 (3,1)+i/3 

(d) 


2-iii-a {de){u){d){ue) 

20 

(3, 2)i/g (1,2),3/2 (3, 1 )+i/3 

(c) 


2-iii-a {de){u){d){ue) 

20 

(3, 2)_i/g (8, 2 )_,_i/2 (3, 1)+i/3 

(d) 


2-iii-b {de){d){u){ue) 

19 

(3, 2)-i/6 (3, 2)+i/6 (3, 1)+i/3 

(d) 


2-iii-b {de){d){u){ue) 

19 

(3, 2)-i/6 (6,2)+i/6 (3, 1)+i/3 

(d) 


2-iii-b {de){d){u){ue) 

20 

(3,2)_3/g (3,1 )_i/3 (3,1)+i/3 

(d) 


2-iii-b {de){d){u){ue) 

20 

(3,2)_3/g (6,1 )_i/3 (3,1)+i/3 

(d) 


4-i {de){u){u){de) 

20 

(3,2)7/g (1,2)3/2 (3,2)+3/g 

(c) 


4-ii-a {uu){d){e){de) 

20 

(6, 1)+4/3 (3, 2)+7/6 (3, 2)+1/6 

(a) 


5-ii-b {ue){e){u){dd) 

19 

(3, l)-l/3 (3, 2)5/g (6, l)2/3 

(b) 

442-i^' 


d 

O 

1—1 

1—1 

BL# 

S S' S" Diagram 

Add. Int. 

2 {ud){ue){de) 

19, 20 

(1,2)+3/2 (3,1)_i/3 (3,2)_i/6 

(d) 


2 {ud){ue){de) 

19, 20 

(8)2)_,_3/2 (3, 1)_i/ 3 (3,2)_3/g 

(d) 



TABLE V: Decompositions that generate d = 7 neutrino mass operator LLHHHH^ at the 2 loop 
level. The topologies of the neutrino mass diagrams in the column of “Diagram” are shown in 
Fia. I j($l For the decompositions in this table, unless some severe fine-tuning of parameters is done, 
the short-range contributions will dominate over the mass mechanism of double beta decay. 
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T-I # Op. 

BL # 

S S' 

2-ii-b {ud){e){u){de) 19, 20 

(8, 2),i/2 (8, 2)i/ 2 (3, 2),i/g 

4-i {de){u){u){de) 

20 

(3,2)-7/6 (8, 2)-i/ 2 (3, 2),i/g 

4-ii-a {uu){d){e){de) 

20 

(6, l)+4/3 (3, l)+5/3 (3, 2),7/g 

4-ii-b {uu){e){d){de) 

20 

(6, l)-H4/3 (6, l)+i/3 (3, 2)+i/g 

4-ii-b {uu){e){d){de) 

20 

(6, l)+4/3 (6, 2),5/g (3, 2),7/g 

5-ii-a {ue){u){e){dd) 

19 

(3, l)-l/3 (6, l)+l/3 (6, l)2/3 

5-ii-a {ue){u){e){dd) 

19 

(3,1)i/ 3 (6,2)i/g (6,1)2/3 

5-ii-b {ue){e){u){dd) 

19 

(3, l)-l/3 (3, l)-4/3 (6, l)-2/3 


T-H # Op. 

BL# 

S S' S" 

4 {uu){de){de) 

20 

(6, l)+4/3 (3,2)7/g (3,2)4/g 

5 {ue){ue){dd) 

19 

(3, l)-l/3 (3, l)-l/3 (6, 1)+2/3 


TABLE VI: Decompositions that generate neutrino masses at 3-loop. Some example diagrams 
are given in the main text. For the decompositions in this table, unless some severe fine-tuning 
of parameters is done, the short-range contributions will dominate over the mass mechanism of 
double beta deeay. 
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T-I# 

Op. 

BL # 

S Ip S' 

3-i 

{uu){e){e){dd) 

- 

(6, l)+4/3 (6, l)+l/3 (6, l)-2/3 

3-ii 

{uu){d){d){ee) 
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{dd){u){u){ee) 
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TABLE VII: Decompositions that generate neutrino masses at 4-loop. For the decompositions in 
this table, the short-range contributions will dominate over the mass mechanism of double beta 
decay. 
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